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Abstract 



o 

We analyze the stability of the structure equations of the vacuum in the brane world models, by 
using both the linear (Lyapunov) stability analysis, and the Jacobi stability analysis, the Kosambi- 



Cartan-Chern (KCC) theory. In the brane world models the four dimensional effective Einstein 
equations acquire extra terms, called dark radiation and dark pressure, respectively, which arise 



> 

, from the embedding of the 3-brane in the bulk. Generally, the spherically symmetric vacuum 



solutions of the brane gravitational field equations, have properties quite distinct as compared 
to the standard black hole solutions of general relativity. We close the structure equations by 
assuming a simple linear equation of state for the dark pressure. In this case the vacuum is Jacobi 
stable only for a small range of values of the proportionality constant relating the dark pressure 



o 

oo 
O 

> 

• rH 
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and the dark radiation. The unstable trajectories on the brane behave chaotically, in the sense that 
after a finite radial distance it would be impossible to distinguish the trajectories that were very 
near each other at an initial point. Hence the Jacobi stability analysis offers a powerful method 
for constraining the physical properties of the vacuum on the brane. 
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I. INTRODUCTION 



The idea of embedding our Universe in a higher dimensional space has attracted a con- 
siderable interest recently, due to the proposal by Randall and Sundrum that our four- 
dimensional (4D) spacetime is a three-brane, embedded in a 5D spacetime (the bulk) 
According to the brane world scenario, the physical fields (electromagnetic, Yang-Mills etc.) 
in our 4D Universe are confined to the three brane. Only gravity can freely propagate in 
both the brane and bulk spacetimes, with the gravitational self- couplings not significantly 
modified. Even if the fifth dimension is uncompactified, standard 4D gravity is reproduced 
on the brane. Hence this model allows the presence of large, or even infinite non-compact 
extra dimensions. Our brane is identified to a domain wall in a 5D anti-de Sitter spacetime. 
In the brane world scenario, the fundamental scale of gravity is not the Planck scale, but 
another scale which may be at the TeV level (for a review of the dynamics and geometry of 
brane world models see [2j). Due to the correction terms coming from the extra dimensions, 
significant deviations from the standard Einstein theory occur at very high energies 
Gravity is largely modified at the electro-weak scale of about 1 TeV. The cosmological and 
astrophysical implications of the brane world theories have been extensively investigated in 
the physical literature {4, ^]. 

The static vacuum gravitational field equations on the brane depend on the generally 
unknown Weyl stresses, which can be expressed in terms of two functions, called the dark 
radiation U and the dark pressure P terms {the projections of the Weyl curvature of the 
bulk, generating non-local brane stresses) 

m. 



Several classes of spherically symmetric solutions of the static gravitational field equa- 
tions in the vacuum on the brane have been obtained in j^, n]]. As a possible physical 
application of these solutions the behavior of the angular velocity v tg of the test particles in 
stable circular orbits has been considered 0, [h], 11, 12]. The observed properties of the ro- 
tational galactic curves [3] can be naturally explained in the brane world models, without 



introducing any additional hypothesis. The galaxy is embedded in a modified, spherically 
symmetric geometry, generated by the non-zero contribution of the Weyl tensor from the 
bulk. The extra terms, which can be described in terms of the dark radiation term U and 
the dark pressure term P, act as a "matter" distribution outside the galaxy. The particles 
moving in this geometry feel the gravitational effects of U, which can be expressed in terms 



2 



of an equivalent mass Mu (the dark mass^ 



12j | . The role of the cross-over length scale 



in the possible explanation of the dark-matter phenomenon in the brane world model was 



investigated in 



14 



also considered in 



. Similar interpretations of the dark matter as a bulk effect have been 



15|- 



For standard general relativistic spherical compact objects the exterior space-time is 
described by the Schwarzschild metric. In the five dimensional brane world models, the high 
energy corrections to the energy density, together with the Weyl stresses from bulk gravitons, 
imply that on the brane the exterior metric of a static star is no longer the Schwarzschild 
metric 6]. The presence of the Weyl stresses also means that the matching conditions do 
not have a unique solution on the brane; the knowledge of the five-dimensional Weyl tensor 
is needed as a minimum condition for uniqueness. 

Static, spherically symmetric exterior vacuum solutions of the brane world models have 
been obtained first in jy] and in The first of these solutions has the mathematical form 
of the Reissner-Nordstrom solution of the standard general relativity, in which a tidal Weyl 
parameter plays the role of the electric charge of the general relativistic solution f|. A 
second exterior solution, which also matches a constant density interior, has been derived 
in Other vacuum solutions of the field equations in the brane world model have been 



obtained in 



13, 



18 



Generally, the vacuum field equations on the brane can be reduced to a system of two 
ordinary differential equations, which describe all the geometric properties of the vacuum 
as functions of the dark pressure and dark radiation terms Q|. In order to close the system 
of vacuum field equations on the brane a functional relation between these two quantities is 
necessary. 

Hence, a first possible approach to the study of the vacuum brane consists in adopting 
an explicit equation of state for the dark pressure as a function of the dark radiation. The 
explicit form of this equation can be obtained by assuming that the brane obeys different 
geometrical or physical conditions. The existence of the group of the homology transforma- 
tions on the brane uniquely fixes the form of the equation of state [8]). By assuming that the 
vacuum admits the group of conformal motions leads to the full determination of the equa- 
tion of state of the dark pressure . By imposing the condition of the constancy of the 
rotational velocity curves for particle in stable orbits closes the system of the field equations, 
and the dark radiation and pressure can be obtained explicitly [8, 
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11. 



12J . This approach 
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has the advantage of making predictions that may be tested directly by observations. 

It is the purpose of the present paper to consider an alternative possibility for constraining 
the equation of state of the dark pressure on the brane. The structure equations of the 
vacuum can be transformed to an autonomous system of two differential equations, which 
in turn may be reduced to a single, second order differential equation. A second order 
differential equation can be investigated in geometric terms by using the general path-space 
;heqry of Kosambi-Cartan-Chern (KCC-theory) inspired by the geometry of a Finsler space 



21 
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23[ . The KCC theory is a differential geometric theory of the variational equations for 
the deviation of the whole trajectory to nearby ones. By associating a non-linear connection 
and a Berwald type connection to the differential system, five geometrical invariants are 



obtained. The second invariant gives the Jacobi stability of the system 



24 



25] . The KCC 



theory has 



sec 



24, 



3een applied for the study of different physical, biochemical or technical systems 
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28 



29| and references therein). 



As a toy model for the applications of the KCC theory for the study of the Jacobi stability 
of the vacuum in the brane world models we consider the case of the linear equation of state 
for the dark pressure, P = jU. The vacuum on the brane is Jacobi stable only for a very 
limited range of the proportionality constant 7, —0.5 < 7 < 0.67, and it is Jacobi unstable 
for all the other possible values of 7. Hence the natural physical requirement of Jacobi 
stability provides a very strong constraint on the equation of state of the dark pressure, and 
on the physical/geometrical properties of the model. Such a constraint cannot be obtained 
from the linear (Lyapunov) stability analysis, which is also considered in detail. 

The present paper is organized as follows. The field equations for the vacuum on the 
brane are written down in Section II. The structure equations of the vacuum are derived in 
Section III. We review the mathematical formalism of the KCC theory in Section IV. The 
linear stability analysis of the structure equations of the vacuum for a linear equation of 
state for the dark pressure is performed in Section V. The Jacobi stability of the structure 
equations is analyzed in Section VI. We discuss and conclude our results in Section VII. 
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II. THE FIELD EQUATIONS FOR STATIC, SPHERICALLY SYMMETRIC VAC- 
UUM BRANES 

In the present Section we briefly describe the basic mathematical formalism of the brane 
world models, present the field equations for a static, spherically symmetric vacuum brane, 
and discuss some of their consequences. 



A. The field equations in the brane world models 

We start by considering a five dimensional (5D) spacetime (the bulk), with a single four- 
dimensional (4D) brane, on which matter is confined. The 4D brane world (^M,g uu ) is 
located at a hypersurface [B (X A ^ = 0) in the 5D bulk spacetime (( 5 'M, <7ab), of which 
coordinates are described by X A ,A = 0, 1, ...,4. The induced 4D coordinates on the brane 
are x^ 1 , /i = 0, 1, 2, 3. 

The action of the system is given by js} 



S — Sbuik + s, 



' brane > 



(1) 



where 



and 



S, 



bulk 



(5)M 



2*1 



(5) i?+ (5) L m + A 5 



S, 



brane 



-(5), 



( 4 )M 



Tf -R' ± + Lbrane (g a /3, 4>) + Af, 



d X, 



(2) 



(3) 



where k\ = 87CG5 is the 5D gravitational constant, and ^L m are the 5D scalar curvature 
and the matter Lagrangian in the bulk, L brane (g a p,ip) is the 4D Lagrangian, which is given 
by a generic functional of the brane metric g a p and of the matter fields ip, is the trace 
of the extrinsic curvature on either side of the brane, and A 5 and (the constant brane 
tension) are the negative vacuum energy densities in the bulk and on the brane, respectively. 
The Einstein field equations in the bulk are given by 



^Gjj = k 2 5 ^T u , ®T„ = -A 5 (5 W + 5(B) [-Xb i5) gu + T„] 



(4) 



where 



(5) T U = -2 



5^g IJ 



(5) 
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is the energy-momentum tensor of bulk matter fields, while T^ v is the energy-momentum 
tensor localized on the brane and which is defined by 

SL brane 

T ^ = - 2 ^r + 3„L brane . (6) 

The delta function 5 (B) denotes the localization of brane contribution. In the 5D space- 
time a brane is a fixed point of the Z 2 symmetry. The basic equations on the brane are 
obtained by projections onto the brane world. The induced 4D metric is gu = ^gij — ninj, 
where rij is the space-like unit vector field normal to the brane hypersurface ^M. In the 
following we assume ®L m = 0. In the brane world models only gravity can probe the extra 
dimensions. 

Assuming a metric of the form ds 2 = (rijnj + gij)dx l dx J , with njdx 1 = d\ the unit 
normal to the \ = constant hypersurfaces and gu the induced metric on \ = constant 
hypersurfaces, the effective 4D gravitational equation on the brane takes the form 

G>i> = ~^g^iv + k^T^ v + k 5 S^ u — E/j, u , (7) 

where S^ u is the local quadratic energy-momentum correction 

11 1 

S»v = Y2 TTflv ~ -^T^Tva + (3T a/3 T Q/3 - T 2 ) , (8) 

and E^ u is the non-local effect from the free bulk gravitational field, the transmitted pro- 
jection of the bulk Weyl tensor Cjajb, Eij = CjAjB nA n B , with the property Ejj — > 
E^SjSj as x — > 0. We have also denoted k\ = 8ttG, with G the usual 4D gravitational 
constant. 

The 4D cosmological constant, A, and the 4D coupling constant, k^ are related by 
A = fc|(A 5 + k 2 \l/6)/2 and k\ = /c|Ab/6, respectively. In the limit A^ 1 -> we recover 
standard general relativity 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation 
of the energy- momentum tensor of the matter on the brane, E V T^ = 0, where D v denotes 
the brane covariant derivative. Moreover, from the contracted Bianchi identities on the 
brane it follows that the projected Weyl tensor obeys the constraint D V E^ = k\D v S^ '. 

The symmetry properties of E^ v imply that_in general we can decompose it irreducibly 
with respect to a chosen 4-velocity field as 



Efja, — —k 



U Uf,u v + -h^ + P^v + 2Q( fi u v ) 



(9) 
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where k = k 5 /k 4 , h^ v = + u^u u projects orthogonal to w M , the "dark radiation" 
term U = —k~ A E^ u u >l u y is a scalar, = k~ 4 h^E a/ 3U^ is a spatial vector and P^ u = 
—k~ 4 \h(n a h v ) 13 — g/i^/i^l E a p is a spatial, symmetric and trace-free tensor. 

In the case of the vacuum state we have p = p = 0, T^ u = 0, and consequently = 0. 
Therefore the field equation describing a static brane takes the form 

R^ = -E^ + A 9fiv , (10) 

with the trace R of the Ricci tensor satisfying the condition R = = 4A. 

In the vacuum case E^ satisfies the constraint D V E^ = 0. In an inertial frame at any 
point on the brane we have w M = 5g and /j™ = diag(0, 1, 1, 1). In a static vacuum = 



and the constraint for E^ v takes the form 



7] 



l -DJJ + hjA^ + D'P^ + A V P^ = 0, (11) 

where = v^DuU^ is the 4-acceleration. In the static spherically symmetric case we 
may chose A^ = A{r)r^ and P^ v = P{r) (r^r u — hh^j, where A{r) and P{r) (the "dark 
pressure") are some scalar functions of the radial distance r, and r M is a unit radial vector 

B. The gravitational field equations for a static spherically symmetric brane 

In the following we will restrict our study to the static and spherically symmetric metric 
given by 

rfs 2 = _ e u{r) dt 2 + e X(r) dr 2 + r 2 ^2 + ^2 ^ 

With the metric given by ffl2|) the gravitational field equations and the effective energy- 



momentum tensor conservation equation in the vacuum take the form 



aa 



e- X (\--)+\ = 3aU + A, (13) 



^"2 iy J ^2 



+ \) -\ = a(U + 2P)-A, (14) 



iy ,y2 / y^-i2 



y^ + ^ + ^-^y a{u - P) . A , (15) 

If 4-2 P' 6P 

(16) 



2U + P r(2U + Py 
where ' = d/dr, and we have denoted a = 16nG/k 4 \b- 



The field equations (fTBl - (TT6]) can be interpreted as describing an isotropic "matter distri- 
bution", with the effective energy density p eS , radial pressure P eff and orthogonal pressure 
Pf, respectively, so that p eS = 3aU + A, P cff = aU + 2aP - A and Pf = aU - aP - A, 
respectively, which gives the condition p eff — P eff — 2Pf l = 4A = constant. This is expected 
for the 'radiation' like source, for which the projection of the bulk Weyl tensor is trace-less, 
J5£ = 0. 

III. STRUCTURE EQUATIONS OF THE VACUUM IN THE BRANE WORLD 
MODELS 

Eq. ffTB"]) can immediately be integrated to give 

— 1 — — — IMnJll _ V, (17) 
r r 3 

where C\ is an arbitrary constant of integration, and we denoted 

GMu(r) = 3a [ U{r)r 2 dr. (18) 
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The function My is the gravitational mass corresponding to the dark radiation term 
(the dark mass). For U — the metric coefficient given by Eq. ( fTTj) must tend to the 
standard general relativistic Schwarzschild metric coefficient, which gives C\ = 2GM, where 
M = constant is the baryonic (usual) mass of the gravitating system. 

By substituting v' given by Eq. (|T6|) into Eq. ([2]) and with the use of Eq. (fT7|) we obtain 
the following system of differential equations satisfied by the dark radiation term U, the 
dark pressure P and the dark mass Mjj, describing the vacuum gravitational field, exterior 
to a massive body, in the brane world model 8]: 

^ = ^W. (19) 
dr G 

dU (2U + P) [2GM + GMp - |Ar 3 + a(U + 2P) r 3 ] dP QP 

dr ~ r 2( 1 _2GM_M 1L _ A r 2) " 2 dr r ' ^ 20 ^ 

V r r 3 / 

To close the system a supplementary functional relation between one of the unknowns U, 
P and Mu is needed. Generally, this equation of state is given in the form P = P(U). Once 
this relation is known, Eqs. (TT9|) (T20T) give a full description of the geometrical properties of 
the vacuum on the brane. 
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In the following we will restrict our analysis to the case A = 0. Then the system of 

equations ffT9l) and fl20l) can be transformed to an autonomous system of differential equations 

by means of the transformations 

2GM CM, ., , . 

q = + n = 3ar 2 U, (21) 

p = 3ar 2 P, 9 = lnr, (22) 

where fi and p are the "reduced" dark radiation and pressure, respectively. 

With the use of the new variables given by Eqs. ( 1211) and ( 1221) . Eqs. ( |T9l) and ( 120]) become 

§ = "-«• < 23 > 

^ = _(2^[ i±i (^ ±M _ 2 | + 
at' 1 — g at* 

Eqs. ( fl~9l) and (l20l) . or, equivalently, ( J23l) and ( |24l) . are called the structure equations of 

the vacuum on the brane [8(. In order to close the system of equations ( 123]) and ( 1241) an 

"equation of state" p = p(/x), relating the reduced dark radiation and the dark pressure 

terms, is needed. 

The structure equations of the vacuum on the brane can be solved exactly in two cases, 
corresponding to some simple equations of state of the dark pressure. In the first case we 
impose the equation of state 2/i + p = 0. From Eq. f l24l) we immediately obtain fi = Qe~ 29 , 
while Eq. ( |23l) gives q (8) = — Qe~ 2e -\-Uoe~ , where Q and Uq = 2GM are arbitrary constants 
of integration. Therefore we obtain the vacuum brane solution 

U = -? = ^\, (25) 
2 3a r 4 ' v ; 



e 



-A 



2GM Q 



+ ^- (26) 



This solution was first obtained in [6], and therefore corresponds to an equation of state 
of the dark pressure of the form P = —2U. The second case in which the vacuum structure 
equations can be integrated exactly corresponds to the equation of state fi + 2p = 0. Then 
Eq. (|24l) gives q^= 2/3, and the corresponding solution of the gravitational field equations 
on the brane is 



U = -2P= — — , e v = C r 2 , e~ x = -. (27) 



1 

, e" = Uqt", e " = - 
9ar 2 3 

This solution corresponds to an equation of state of the dark pressure of the form P 
-U/2. 



IV. KOSAMBI-CARTAN-CHERN (KCC) THEORY AND JACOBI STABILITY 



We recall the basics of KCC-theory to be used in the sequel. Our exposition follows [24]. 

Let M. be a real, smooth n-dimensional manifold and let TM. be its tangent bundle. Let 
(x l ) = (x 1 , x 2 , x n ), (y l ) = (y 1 , y 2 , y n ) and the time t be a 2n + 1 coordinates system 
of an open connected subset Q of the Euclidian (2n + 1) dimensional space R n x R n x R 1 , 
where 

/ dx 1 dx 2 dx n \ , . 

(*-■*-•■--*)■ (28) 

We assume that t is an absolute invariant, and therefore the only admissible change of 
coordinates will be 

t = t, 5? = x* (x\x 2 , ...,x n ),ie {1, 2, n} . (29) 

The equations of motion of a dynamical system can be derived from a Lagrangian L via 
the Euler-Lagrange equations, 

d dL dL , . 

*ajji-*f =F,,, = 1 *-"- nt (30) 

where Fi, i = 1, 2, n, is the external force 27|. The triplet (M, L, Fj) is called a Finslerian 



mechanical system 



261 ] . For a regular Lagrangian L, the Euler-Lagrange equations given by 



Eq. ( pOll are equivalent to a system of second-order differential equations 

" ' ; +2G l (x\yi,t)=0,ie{l,2,...,n}, (31) 



where each function G l (x\y : ',t) is C°° in a neighborhood of some initial conditions 
((x)o j (y)o > *o) m ^- The system given by Eq. (l3"Ij) is equivalent to a vector field (semispray) 
S, where 

5 = ^- 2G W,*)|l> (32) 



which determines a non-linear connection iVj defined as 



27] 



iVj = (33) 
3 dyi 

More generally, one can start from an arbitrary system of second-order differential equa- 
tions on the form (131j) . where no a priori given Lagrangean function is assumed, and study 
the behavior of its trajectories by analogy with the trajectories of the Euler-Lagrange equa- 
tions. 
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For a non-singular coordinate transformations given by Eq. (|29l) . we define the KCC- 
covariant differential of a vector field £ l (x) on the open subset Q C R n x R n x R 1 as 



24 



25 



2g,|2S| 



r>e = de 

dt dt 



For P = y l we obtain 



Dy i 



My j - 2G l 



(34) 
(35) 



dt ' 3 

where the contravariant vector field e l on Q is called the first KCC invariant. 

Let us now vary the trajectories x l (t) of the system (13T|) into nearby ones according to 



x l (t) =x i (t)+r)C{t), 



(36) 



where \tj\ is a small parameter and £*(£) are the components of some contravariant vector 
field defined along the path x l {t). Substituting E q s . f|36l) into Eqs. ( 13TI) and taking the limit 



r] — > we obtain the variational equations 



24 



25 



28 



29] 



d 2 ^ dP 3G { 

ZJL + 2N l — + 2—P 
dt 2 + > dt +Z dx^ 



0. 



(37) 



By using the KCC-covariant differential we can write Eq. fl37|) in the covariant form 

d 2 c 



dt 2 



where we have denoted 



P\ -2^ -2G l Gl, + y l ^ + NlN l j + 1 



^ 2GW +/ N] 
and G % ai = dNUdy 1 is called the Berwald connection 



dt 



(38) 



(39) 



24 



25 



23, |28j. Eq. flU is called 



the Jacobi equation, and Pj is called the second KCC-invariant or the deviation curvature 
tensor. When the system ( |3T1) describes the geodesic equations in either Riemann or Finsler 
geometry, Eq. ( |38j) is the Jacobi field equation. 



The third, fourth and fifth invariants of the system fl3T|) are given by 29] 



P l = - 



1 (dP\ dPl 



pi 



dP, 



dG 



.1^ 



(40) 



Qyk Qyj J ' J kl Qyl ' 3 k l Qyl 

The third invariant is interpreted as a torsion tensor, while the fourth and fifth invariants 
are the Riemann- Christoffel curvature tensor, and the Douglas tensor, respectively j^J. In 
a Berwald space these tensors always exist, and they describe the geometrical properties of 
a system of second-order differential equations. 
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In many physical applications we are interested in the behavior of the trajectories of the 
system (1311) in a vicinity of a point x l (to), where for simplicity one can take to = 0. We 
will consider the trajectories x l = x l (t) as curves in the Euclidean space (R n , (•,•)), where 
(., .) is the canonical inner product of R n . As for the deviation vector £ we assume that it 
satisfies the initial conditions £ (0) = O and £ (0) = W ^ O, where O G R n is the null vector 



24 



25|. 



For any two vectors X, Y e R n we define an adapted inner product ((.,.)) to the deviation 
tensor £ by ((X,Y)) : = 1/ (W, W) ■ (X,Y). We also have \ \W\\ 2 := ((W,W)) = 1. 

Thus, the focusing tendency of the trajectories around to = can be described as follows: 
if ||£(t)|| < t 2 , t ps + , the trajectories are bunching together, while if ||£(t)|| > t 2 , t ~ + , 



the trajectories are dispersing 



24 



251 ] . In terms of the deviation curvature tensor the 



focusing tendency of the trajectories can be described as follows: The trajectories of the 
system of equations (|3~T!) are bunching together for t ~ + if and only if the real part of the 
eigenvalues of Pj (0) are strictly negative, and they are dispersing if and only if the real part 
of the eigenvalues of Pj (0) are strict positive [2J, l25j ]. 

Based on these considerations we can define the Jacobi stability for a dynamical system 



as follows 



25 



29|: 



Definition: If the system of differential equations (}3~Tj) satisfies the initial conditions 
\\ xl (*o) ~ x% (^o)ll — 0, \\x l (t ) — x % (to) 1 1 7^ 0, with respect to the norm ||.|| induced by a 
positive definite inner product, then the trajectories of (I3T1) are Jacobi stable if and only if 
the real parts of the eigenvalues of the deviation tensor Pj are strictly negative everywhere, 
and Jacobi unstable, otherwise. 

The focussing behavior of the trajectories near the origin is represented in Fig. [TJ 



V. LINEAR STABILITY ANALYSIS OF THE VACUUM STRUCTURE EQUA- 
TIONS ON THE BRANE FOR THE LINEAR EQUATION OF STATE OF THE 
DARK PRESSURE 



Since generally the structure equations of the vacuum on the brane cannot be solved 
exactly, in this Section we shall analyze them by using methods from the qualitative analysis 



of dynamical systems 111 ], by closely following the approach of 



30l . We consider the case 



in which the dark pressure is proportional to the dark radiation, P = jU, where 7 is an 
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2 <t\ t^0+ 



|£(t)|| 2 >t 2 , t*0+ 



FIG. 1: Behavior of trajectories near zero. 



arbitrary constant, which can take both positive and negative values. As we have seen in 
Section 1111} several classes of exact solutions of the vacuum gravitational field equations on 
the brane can be described by an equation of state of this form. In the reduced variables /i 
and p the linear equation of state is p = 7/i, and the structure equations of the gravitational 
field on the brane have the form 



dq 

:i+27) ^ 2(1 - 7) ,- (^M«+^] 

do 1 — q 



(41) 



(42) 



Let us firstly analyze the special case where 7 = —1/2. Then, by virtue of the second 
Eq. ( )42l) . we obtain two possible exact solutions, either q = {i = or q = fi = 2/3. The 
first of these solutions correspond to the vanishing of all physical quantities, and therefore 
we can discard it as unphysical. The second case corresponds to an exact solution, which 
has been discussed in Section IIHI 

Let us henceforth assume that 7 7^ —1/2, and rewrite the system of equation into the 
following form, 

(43) 



dq 

de = - q + ^ 



dfi = 2(1 - 7) 7 + 2 fi [q + gg>] 
d6 1 + 2 7 ^ 1 + 27 1 - q 



which is finally written as 



dd 



(44) 



(45) 
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where we have denoted 



r 


|. A-( 







2(l- 7 )/(l + 2 7 ) 





\l+2 7 " 1-q 



7+2 /^[g+^^Mj 



(46) 



The system of equations (1451) has two critical points, X = (0, 0), and 



3(1-7) 3(1-7) 
v 7 2 + 7 + 7' 7 2 + 7 + 7 / 

For 7 = 1, the two critical points of the system coincide. Depending on the values of 7, 
these points lie in different regions of the phase space plane (q,fj). 

Since the term ||-B||/||£|| — ^ as 1 1^| | — > 0, the system of equations f|4"5]) can be linearized 
at the critical point X . The two eigenvalues of the matrix A are given by 77 = —1 and 
T2 = 2(1 — 7 )/(l + 27), and determine the characteristics of the critical point Xq. For 
7 G (—00,— 1/2) U [l,oo) both eigenvalues are negative and unequal. Therefore, for such 
values of 7 the point X is an improper asymptotically stable node. 

If 7 G (—1/2,1), we find one positive and one negative eigenvalue, which corresponds 
to an unstable saddle point at the point X . Moreover, since the matrix dA/d!;(X Q ) has 
real non-vanishing eigenvalues, the point X is hyperbolic. This implies that the properties 
of the linearized system are also valid for the full non-linear system near the point Xq. It 
should be mentioned however, that this first critical point is the less interesting one from a 
physical point of view, since it corresponds to the 'trivial' case where both physical variables 
vanish. 

As we have seen in Section IHIt the structure equations can also be solved exactly for the 
value 7 = —2. In that case, the non-linear term B in Eq. (f45|) identically vanishes, and the 
system of equations becomes a simple linear system of differential equations. For 7 = —2 
the two eigenvalues of A are given by 77 = —1 and r 2 = —2, and the two corresponding 
eigenvectors are linearly independent. The general solution can be written as follows 



e 7= _ 2 = (g + /io) \e- e + /i e" 2e , (48) 




where go — g(0) and /io = /J-(O) . One can easily transform this solution back into the radial 
coordinate r form by using 9 = ln(r), thus obtaining 

/i 7 =- 2 = ^, g 7 =- 2 = - + //ol --^)- (49) 
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7 


— oo 


-0.5 0.67 1 +oo 


r± 


real 


complex real real 




Re r+ < 


r + 


+ - - 


r _ 


- | + 




Saddle 


Stable spiral Stable node | Saddle 



TABLE I: Linear stability of the stable point Xy. 

Let us now analyze the qualitative behavior of the second critical point X 1 . To do this, 
one has to Taylor expand the right-hand sides of Eqs. (1431) and (jHJ) around X 7 and obtain the 
matrix A which corresponds to the system, linearized around X 7 . This linearization is again 
allowed since the resulting non-linear term, h say, also satisfies the condition ||^||/||£|| — > 
as ll£ll — > -^7- The resulting matrix reads 



A 



-1 1 

3(-7 2 +57-4) 7-1 
,(2+7)2(1+27) 2+ 7/ 



(50) 



and its two eigenvalues are given by 

-3 - 67 ± v/167 4 + 87 s + 132 7 2 - 287 - 47 



47 2 + IO7 + 4 



(51) 



For —0.5 < 7 < 0.674865 the argument of the square root becomes negative and the 
eigenvalues complex. Moreover, the values 7 = — 1/2 and 7 = —2 have to be excluded, since 
the eigenvalues in Eq. (1501) are not defined in these cases. However, both cases have been 
treated separately above. 

If 7 G (—00, — 1/2) U (1, 00), then X 1 corresponds to an unstable saddle point and for 7 G 
(0.674865, 1) it corresponds to an asymptotically stable improper node. More interesting is 
the parameter range 7 G (—1/2, 0.674865), where the eigenvalues become complex, however, 
their real parts are negative definite. Hence, for those values we find an asymptotically stable 
spiral point at X T Since this point is also a hyperbolic point, the described properties are 
also valid for the non-linear system near that point. 

The behavior of the trajectories is shown, for 7 = — 1 and 7 = 0.4, in Fig. [21 The figures 
show the attracting or repelling character of the steady states, respectively. The results of 
the linear stability analysis of the critical points X 7 are summarized in Table [H 
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1 \ w 



0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 



FIG. 2: Behavior of the trajectories of the structure equations of the vacuum on the brane near 
the critical points for 7 = — 1 (left figure) and 7 = 0.4 (right figure). 

VI. JACOBI STABILITY ANALYSIS OF THE VACUUM GRAVITATIONAL 
FIELD EQUATIONS IN THE BRANE WORLD MODELS 

Since from Eq. (jUJ) we can express fj, as [A — q+ dq/d9, upon substitution in Eq. f|42|) 
we obtain for q the following second order differential equation 



d 2 q 



+ 



d6 2 3 (1 + 2 7 ) (1 - q) 



6 (7 - 1 ) q + 2 ( 7 2 + 7 + 7) q 2 + 3 (4 7 - 1) ^ 



+ (4 7 2 + 7 + 13) g| + (2 7 2 + 5 7 + 2) g 



(52) 



0. 



which can now be studied by means of KCC theory. 

By denoting x = q and dq/dd = dx/d6 = y, Eq. ( 1521) can be written as 

d 2 x 



d6 2 



+ 2G 1 (x,y)=0, 



(53) 



where 



G 1 (x,y) 



6 (7 - 1) x + 2 ( 7 2 + 7 + 7) x 2 + 3 (4 7 - 1) y 



(54) 



6 (1 + 2 7 ) (1 - x) 
+ (4 7 2 + 7 + 13) xy + (2 7 2 + 5 7 + 2) y 2 

As a first step in the KCC stability analysis of the vacuum field equations on the brane 
we obtain the nonlinear connection N± associated to Eq. ( 1521) . and which is given by 
dG 1 3 (4 7 - 1) + (4 7 2 + 7 + 13) x + 2 (2 7 2 + 5 7 + 2) y 



dy 



6 (1 + 2 T ) (1 - x) 



(55) 
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The Berwald connection can be obtained as 

ON} 2 7 2 + 57 + 2 



G\ 



11 dy " 3(1 + 27) (1-x)' 



Finally, the second KCC invariant or the deviation curvature tensor defined as 



(56) 



(57) 



reads now 



Pi (x, y) 



27-2 [61 + 577 + 4 7 2 (9 + 27)] x + 3 (5 + 7 ) 2 x 2 - 2 (2 + 7) (1 + 2 7 ) (5 + 4 7 ) y 



12(l-x) 2 (l + 2 7 ) : 



(58) 



Taking into account that x = q and y — jj, — q, we obtain Pi in the initial variables as 



27-6 (17 + 87 + 2 7 2 ) g + 3 (5 + 7 ) 2 g 2 - 2 (2 + 7) (1 + 2 7 ) (5 + 4 7 ) /i 



12(l-g) 2 (l + 2 7 ) 2 
Evaluating Pi (q, /x) at the critical point X 7 , given by Eq. (HTj) . ew obtain 

4(2 + 7 ) 2 (l + 2 7 ) 

The plot of the function Pi (X 7 ) is represented in Fig. [31 

Pi 1 

(40 ; 

; 
10 ; 
10 1 



-10 



-6 




r 



(59) 



(60) 



FIG. 3: The deviation curvature tensor Pi (X~) as a function of 7. 

Using the discussion of Section |V] our main results on the Linear stability and Jacobi 
stability of the critical point X 7 of the vacuum field equations in the brane world models 
can be summarized in Table [Til 
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7 


— oo 


-0.5 0.67 


1 


+oo 




+ 


1 " 1 


+ 


+ 


Linear stability 


Saddle 


Stable 


Stable 


Saddle 


of X 7 


point 


spiral 


node 


point 


Jacobi stability 


Jacobi 


Jacobi 


Jacobi 


Jacobi 


of X 7 


unstable 


stable 


unstable | 


unstable 



TABLE II: Linear and Jacobi stability of the stable point 



VII. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the stability properties of the vacuum gravita- 
tional field equations in the brane world models. For the analysis of the stability we have 
used two methods, the Lyapunov (linear) stability analysis, and the so-called Jacobi stabil- 
ity analysis, or the KCC theory. The study of the stability has been done by analyzing the 
behavior of the steady state X 1 of the structure equations of the vacuum on the brane. The 
Lyapunov stability analysis involves the linearization of the dynamical system via the Jaco- 
bian matrix of a non-linear system, while the KCC theory addresses the Lyapunov stability 
of a whole trajectory in a tubular region 24]. 

By using the KCC theory we have shown that the vacuum on the brane is Jacobi unstable 
for most of the values of the parameter 7. The stability region is reduced to a very narrow 
range of 7, 7 G (—0.5, 0.67). For all other values of 7 the vacuum on the brane is unstable, 
in the sense that the trajectories of the structure equations will disperse when approaching 
the origin of the coordinate system. 

In a previous paper ({24]) we have regarded the Jacobi stability of a dynamical system 
as the robustness of the system to small perturbations of the whole trajectory. This is 
a very convenable way of regarding the resistance of limit cycles to small perturbation of 
trajectories. On the other hand, we may regard the Jacobi stability for other types of 
dynamical systems (like the one in the present paper) as the resistance of a whole trajectory 
to the onset of chaos due to small perturbations of the whole trajectory. This interpretation is 
based on the generally accepted definition of chaos, namely a compact manifold M on which 
the geodesic trajectories deviate exponentially fast. This is obviously related to the curvature 
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of the base manifold (see section IV). The Jacobi (in)stability is a natural generalization 
of the (in) stability of the geodesic flow on a differentiable manifold endowed with a metric 
(Riemannian or Finslerian) to the non-metric setting. In other words, we may say that 
Jacobi unstable trajectories of a dynamical system behave chaotically in the sense that after 
a finite interval of time it would be impossible to distinguish the trajectories that were very 
near each other at an initial moment. 

We have found that there is a good correlation between the linear stability of the critical 
point X 7 and the robustness of the corresponding trajectory to small perturbations. Indeed, 
for small values of the parameter 7 the saddle point is also Jacobi unstable (7 < —0.5 and 
7 > 1), while the stable spiral obtained for —0.5 < 7 < 0.674865 is also robust to a small 
perturbation of all trajectories. 

It is interesting to remark that for the interval 0.674865 < 7 < 1, the stable node is 
actually Jacobi unstable. In other words, even the system trajectories are attracted by the 
critical point X 7 one has to be aware of the fact that they are not stable to small perturbation 
of the whole trajectory. This means that one might witness chaotic behavior of the system 
trajectories before they enters in a neighborhood of X 1 . We have here a sort of stability 
artifact that cannot be found without using the powerful method of Jacobi stability analysis. 
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